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00 ' A superglass is a phase of matter which is characterized at the same time by superfluidity and a 

f^ ' frozen amorphous structure. We introduce a model of interacting bosons in three dimensions that 

("^ , displays this phase unambiguously and that can be analyzed exactly or using controlled approxima- 

04 ' tions. Employing a mapping between quantum Hamiltonians and classical Fokker-Planck operators, 

we show that the ground state wavefunction of the quantum model is proportional to the Boltzmann 
O , measure of classical hard spheres. This connection allows us to obtain quantitative results on static 

and dynamic quantum correlation functions. In particular, by translating known results on the 
glassy dynamics of Brownian hard spheres we work out the properties of the superglass phase and 
of the quantum phase transition between the superfluid and the superglass phase. 
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P^ ; I. INTRODUCTION 

C/3 . Solids do not flow. This apparently tautological statement can be wrong in two ways. First, solidity is in general a 

timescale-dependent phenomenon. Crystal or glasses, well known solids on most experimentally relevant timescales, 

do flow if one waits long enough (see Ref. 1 for example). Second, at very small temperatures, quantum solids may 

d • become superfluids as suggested theoretically in the early seventies^"'*. A striking property of this phase of matter is 

the non-classical rotational inertia: a supersolid placed at thermal equilibrium in a rotating container does not rotate, 

I ' at least not completely, for small angular velocity lo. In particular, its angular momentum is reduced from its classical 

' ^ ' value /ci w by a fraction fs which is called superfluid fraction. 

^ The first promising experimental evidence of supersolidity was found only four years ago by Kim and Chan^'^ 

through measurements of a non-classical rotational momentum of inertia (NCRI) in solid He*^ at low temperature. 
Unfortunately, the physical mechanism underlying the observed effect is still unclear. In particular, despite intense 
experimental activity, it has not been established that the observed NCRI in solid He^ is accompanied by superflow 
►-^ " (see, however, the very recent work reported in Ref. 7). 

(^ ■ Nonetheless, some agreement has been reached by now on few basic facts*: theoretical works have revealed that 

ly-s I equilibrium He** crystals at zero temperature are commensurate and not supersolids^"^^, or at least have an extremely 
■^ . small superfluid fraction, much too small to account for experimental findings (see also Ref. 12). Experiments 
O^ ' have shown that the supersolid critical temperature and the superfluid density depend sensitively on the detailed 
r — [ , preparation history of the solid samples^^ and on the presence of minute fractions of He'^ impurities^'^. The more 
(•^ ' disordered are the He'' solids the larger is the supersolid signal. It can reach surprisingly large values when the solid 
00 ! is formed by a rapid freezing from the normal phase. It can be reduced to an unobservable level by annealing^'*. It 
^D is by now very plausible that rapid freezing produces highly disordered samples and that disorder is at the origin 
^ of the experimental flndings and of "supersolid" signals. Different origins have been highlighted. A recent work 
in Ref. 15 suggested that the core of dislocations could be superfluid. A phenomenological theory of the role of 
dislocations can be found in Ref. 16. Following up Ref. 17 a scenario to explain the experimental findings in terms of 
dislocation induced superfluidity has been put forward in Ref. 18. From the experimental point of view, the role of 
grain boundaries has been analyzed in Ref. 19. Another idea introduced in Refs. 20,21 is that rapid freezing produces 
amorphous glassy solids that become supersolid at low temperature. This "superglass" phase has been detected in 
Path-Integral Monte-Carlo (PIMC) numerical simulations^". These results are encouraging but addressing real-time 
dynamics issues such as the stability and the dynamics of this phase is out of reach of PIMC. Therefore, it is desirable 
to have a complementary analytical approach that can give direct insight into these problems: this is precisely the 
purpose of this work. 

Our strategy will be to focus on a model that displays a superglass phase and that can be concretely analyzed in some 
instances exactly and in others using accurate and controlled approximations at zero temperature. The model consists 
of bosonic particles interacting via a short-range potential in three dimensions. Our approach is based on a mapping 
between quantum Hamiltonians and (classical) Fokker-Planck operators which allows us to obtain results on ground 
state properties and time dependent correlation functions from the analysis of the stochastic dynamics of a classical 
equilibrium system. This connection, already well-known for a few decades^^'^^, has been quite used and further 
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studied in the context of the so caUed Rokshar-Kivelson points^^"^^. We wih mainly focus on systems of interacting 
bosons whose classical counterpart (via the above mapping) is a system of hard spheres. These are known to undergo 
a crystal^^ or a glass^^"'^^ transition at high density. The glass transition takes place when crystallization is avoided 
(either by fast compression or using different types of particles) . The usefulness of the classical-quantum mapping is 
that it allows us to obtain highly non-trivial, very accurate or, at instances, exact results on the original quantum 
many-body problem by translating all the knowledge on the Brownian dynamics of hard spheres. In particular we 
will be able to work out in detail the properties of the superglass, supercrystal and superfluid phases. The drawback 
of this mapping is that the resulting quantum solids are quite different from real He^ crystals, mainly because of the 
much smaller zero-point motion and of some peculiar properties of the excitation spectrum to be discussed below. 

The results of our work go beyond the physics of solid He'* because the superglass phase can emerge also in other 
physical contexts, in particular strongly interacting bosonic mixtures of cold atoms. Furthermore, another motivation 
for our study is to analyze a quantum glassy phase not induced by quenched disorder. There are in fact very few 
examples^^^^^ of that, especially not relying on mean-field or large- A'^ approximations. 

II. SUMMARY OF THE MAIN RESULTS 

In the following section the mapping between classical Fokker-Planck operators and quantum Hamiltonians is 
presented. In particular we will obtain the relation between the classical potential and its quantum counterpart. We 
will work it out for identical bosons as well as for mixtures of bosons (mixtures may be relevant for cold atomic 
systems). Furthermore, we shall introduce the model that we will focus on in this paper, which is the quantum 
counterpart of classical Brownian hard spheres. The quantum potential is still a hard sphere one but it has in 
addition a sticky part at contact (and also a 3-body interaction at contact). 

In section IV we shall work out the zero temperature phase diagram of the quantum model as a function of the 
density p. Actually we will use instead of the density the packing fraction cj) = tt/jct^/G (where a is the sphere 
diameter) which is a more natural variable for hard spheres. We find three possible phases: superfluid, supercrystal 
and superglass. At low packing fraction the system is superfluid but at <}() ~ 0.49 it has a first order transition toward 
a supercrystal phase. The transition toward the superglass phase takes place at higher packing fractions, ((> ~ 0.62. 
It can be achieved only if the transition toward the crystal is avoided either by considering mixtures or compressing 
fast the liquid phase. We will obtain quantitative results for the condensate fraction. In particular for the superglass 
phase we will introduce and compute an observable, akin to the Edwards-Anderson parameter for spin glasses, that 
quantifies the inhomogeneity of the condensate wavefunction in space. The values of the condensate fraction turn out 
to be very small. This is likely due to the particular model we have chosen which admits as ground state wavefunction 
the Gibbs measure of hard spheres. However, exact results guarantee that, although small, the condensate fraction is 
non-zero, i.e. it is not an artifact. In Section V we focus on the real time quantum dynamics close to the transition 
from the superfluid to the superglass. In particular we obtain the behavior of the density correlations as a function 
of time as well as the time dependent condensate fluctuations. Both correlation functions show a similar behavior. 
Approaching the superglass from the superfluid, one finds that after a rapid oscillating decay, there is a plateau 
indicating frozen amorphous density and condensate fluctuations. On much larger timescales these frozen fluctuations 
relax. The transition to the superglass phase takes place when the timescale for this second relaxation goes to infinity 
(or, practically, becomes larger than any relevant experimental timescale.). This behavior is actually very reminiscent 
of the one of classical supercooled liquids approaching their structural glass transition. In Section VI we will generalize 
our model in order to make it more realistic. In fact, within the classical-quantum mapping approach, the energy of 
the ground state is always zero. In order to discuss the behavior of the pressure, the order of the transitions and the 
superfluid properties we shall add a small perturbation to the quantum potential originally introduced in Section III. 
This will allow us to show that the transition from the superfluid to the superglass is unusual to the extent that it 
is thermodynamically first order but has also some characteristics of second order transitions. Finally, conclusions, 
perspectives and relations to experiments and numerical simulations will be presented in Section VII. 

III. FROM CLASSICAL LANGEVIN PARTICLES TO ZERO TEMPERATURE INTERACTING 

BOSONS 

The superglass is an amorphous quantum many-body state of interacting bosons. To find such a state starting from 
a generic quantum Hamiltonian of interacting particles is a daunting task. In this paper, we take an approach that 
allows us to argue for such a phase in a well controlled way, by constructing a local quantum many-body Hamiltonian 
whose ground state is known exactly and can be argued to be a glassy state at large enough densities. This formulation 
is motivated by a generic result stating that for systems without a sign problem there is a simple connection between 



quantum Hamiltonians and the stochastic dynamics of a classical system. This can be used to construct interesting 
Hamiltonians with real and non-negative ground state wavefunctions related to classical equilibrium Boltzmann-Gibbs 
measures. The connection is far-reaching since it allows one to obtain controlled and highly non-trivial results on 
the phase diagram and dynamical properties of a quantum many body problem^^. Moreover, it has been used to 
construct an efhcient algorithm for Quantum Monte Carlo at zero temperature^^. 

Here we shall explore this classical-quantum connection for the particular case of bosonic point particles, following 
the standard route for mapping classical Langevin dynamics for a many-particle system and its associated Fokker- 
Planck operator to a Schrodinger operator^^'^^. 

A. Langevin dynamics: the Fokker-Planck and Schrodinger operators 

Consider N particles whose evolution is determined by the following Langevin equations: 



l^-^ =-^C/Ar(xi,...,XAr)+rj,(i) 



2-1, ...,7V 



(1) 



where ji are friction coefficients, Tjfit) is a Gaussian white thermal noise with variance {vt(t)'>lj {t')) = 
2Tji 6ij Sap 5(t — t'). Furthermore, T is the temperature (with fcs = 1) and a and (3 run from 1 to the spatial 
dimension d (henceforth the boldface notation indicates vectors). The potential will eventually be assumed to be the 
sum of (symmetric) pair potentials. 



C/Ar({x}) = C/Ar(xi, . . . , Xat) = X X! ^«J (^' ^ ^j) 



(2) 



i¥=0 



with Vij = Vji, and g|-C/Ar = V^C/at = Ej(^i) ^Viji^i - ^j)- 

It is well known^^'^^ that the evolution equation for the probability distribution P({x}) can be written as a 
Schrodinger equation in imaginary time: 



dtP 



-HppP 



(3) 



where the Fokker-Planck operator reads: 



^ 1% 9xj 



dx. 



(4) 



The Fokker-Planck operator is non-Hermitian and can be proven to have all eigenvalues larger than or equal to 
zero^^'^"^. The zero eigenvalue corresponds -as it can be readily checked- to the stationary distribution which is the 
equilibrium Gibbs probability measure: 



Pg({x}) = — e'T^^wd^}) = ^ (>-ttY.^^,v,,(::,l,-^j 
Zn Zn 



(5) 



Setting for simplicity h — 1, the Fokker-Planck operator can be mapped into a Hermitian quantum Hamiltonian by 
a similarity transformation 



H ^ e^^'^ Hpp e-^^", 



(6) 



that leads to 
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This expression corresponds to a Hamiltonian for particles with mass nii ~ "/i/(2T) and an effective potential which 
is the sum of a two body and three-body interaction: 



V^({x}) = ^l 



-1^'Un + ^{^^UNf 






(8) 



The eigcnfunctions of the quantum Hamiltonian and of the Fokker Planck operator are in a one to one correspondence. 
Indeed, by applying the similarity transformation introduced above one finds: 



^e{M) « ew Pb({x}) 



(9) 



where Pg indicates the right eigenfunction of the Fokker-Planck operator with eigenvalue E, and '^ e its counterpart 
associated to the quantum Hamiltonian. This also implies that all the eigenvalues E corresponding to the Fokker- 
Planck operator are identical to the ones of the quantum Hamiltonian. 

In particular, this relation, together with Eq. (5), allows one to obtain straightforwardly the ground state wave- 
function of the quantum problem, which is of the Jastrow form^^"'*^ : 



*g({x}) = TT^Iw) 



1 



: exp 



^T.^^3i^^ 



Xj) 



i#J 



(10) 



The logic of the approach we pursue hereafter is the following: we take as starting point Hamiltonians with many- 
body potentials of the form Eq. (8), for which the Jastrow form Eq. (10) is exact. In general, wavefunctions of this 
form lead to more than two-body interactions VAr({x}) (note that also He** has weak higher order interactions). The 
important point is that if the two-body potentials Vij(xi — Xj) are short-ranged {i.e. local) potentials, then Vjv({x}) 
is also local, and thus the many-body Hamiltonians on which we focus are local (non-local Hamiltonian may lead to 
pathological behaviors). 

Because we know exactly the ground state wavefunction, and it is related to a Boltzniann-Gibbs measure for 
a classical system, quantum static correlation functions can be computed in terms of classical static correlation 
functions'*^. Furthermore, as we shall show and noticed by Henley^^, the mapping generalizes also to dynamical 
correlation functions. Hence, we will obtain quantum dynamical correlation functions at zero temperature by analytic 
continuation of classical (stochastic) dynamical correlation functions. 

To simplify the notations, in the following we will fix T = 1 and (7) ~ N^^ ^^ ji = 1. Together with h = 1, this 

fixes the units in both classical and quantum problems. As a consequence the masses in the quantum problem read 

_ 1 jri_ 
m^ - 2 (^) • 

B. Identical bosons 

Let us first consider the simplest case of N identical hosons characterized by the Hamiltonian H in Eq. (7) with 
7i = 7 = 1 and Vij = V . It is important to remark that since the particles are bosons one has to consider only 
many-body states that are completely symmetric under permutation of particles, and study only observables that are 
invariant under permutation of particles (e.g. the density-density correlator). This is clearly not a difficult constraint 
to handle since the Jastrow form (10) with Vij = T^ is completely symmetric. Furthermore, even in the study of 
dynamical correlations (section V) this will not be a problem because if one starts from a probability law completely 
symmetric under permutation of particles, the symmctrization carries over to all later times. This follows trivially 
from the Fokker-Planck evolution Eq. (3), since if the state P and the operator Hpp are both symmetric under 
exchange of particles, so is the time derivative dtP and thus the many-body state thereafter. 

For a given classical isotropic 2-body potential T^(x) = l/(|x|), the resulting quantum potential energy will have 
2-body and 3-body interactions: 

VAr({x}) = Y. ^''"'(^'' - ^^■) + E «'"''°''nx^ - x.> x^ - ^y) 



z;P-^(x) = -VV(x) + \ [Vy(x)]2 



-'L-^Y[r)-V"{r) + \\V\r)f 



(11) 



V 



3-b°dy(x^x') ^ 1 vy(x) . vy(x') = i ^ • ^ r (r) F'(/) , 



where d is the spatial dimension and r = | x| 



C. Binary mixtures of bosons 



It is also worthwhile considering binary mixture of bosons. In fact the superglass phase we shall discuss in the 
following will emerge only when crystallization is avoided. We expect, as it is the case for classical systems, that 



this can be obtained either by fast compression of the hquid state or by considering mixtures of bosons. Mixtures 
typically induce a frustration effect on the crystalline phase whose energy, as a consequence, increases. Instead the 
glass phase, thanks to its disordered structure, is less affected and, hence, starts to compete in terms of energy (or 
dynamical basin of attraction) with the crystalline phase. This conclusion, which we expect to be correct on physical 
grounds, is also suggested by approaches based on variational wavefunctions of the Jastrow form Eq. (10). These 
types of wavefunctions map the quantum problem into a classical one and are very effective in describing a number 
of the observed phenomena in the quantum Bose fluid. Lennard- Jones type potentials (as classical potentials in the 
exponent of the Jastrow wavefunction) have been considered by McMillan^'*, Francis, Chester and Reatto^^ in the 
study of He** superfluids. For a single species of particles, such interaction potentials lead to classical liquid states 
and, at very high densities, might lead to a crystalline state^^. Within Jastrow-like approaches, thanks to this analogy 
with classical systems, we expect the physical effect induced by mixtures to be similar to what happens in classical 
liquids. 

As a consequence, we conjecture that for a binary mixture, a quantum glassy state may emerge at high enough 
densities. Experimentally, the binary mixture can be obtained using two species of bosonic atoms, A and B. It may 
be the case that we do not need to develop a condensate fraction for both species A and B; it would suffice to have a 
superglass of, say, species A, while species B always remains normal, with its role solely that of enabling the formation 
of the glassy state. These type of binary mixture can be obtained in cold atom experiments. 

The generalization of the results of the previous section to a binary mixture is presented in Appendix A. In the 
following we will focus on the identical boson system but one must keep in mind that to stabilize the glassy phase one 
may have to consider mixtures or fast compression of the liquid state. Thus, our strategy will be to study directly 
the glass phase of identical particles and translate the results to the mixture case where this phase actually exists. 

D. Quantum model with hard sphere wavefunction 

Let us investigate what quantum bosonic system in particular corresponds to a classical hard sphere problem upon 
the mapping above. The motivation for this study is that the classical hard sphere packing is a well studied problem in 
the context of glasses for large enough packing fractions^^"^^ , and we will be able to use these results for substantiating 
the notion of the superglass state. We shall focus on a classical potential V{r) = Voexp(— A[(r/CT)^ — 1]), where a is the 
characteristic sphere size, perform the mapping and take the A — > oo limit that enforces the hard sphere constraint. 
Note that although in principle the form of the potential is not very important as long as it is infinite for r < a and 
zero for r > a (in the A — > oo limit), the previous form makes the following discussion particularly easy. We also set 
Vo = 1 and a — I for simplicity: indeed, even if the units have already been fixed, Vq is irrelevant in the A ^ oo limit, 
while a change in a is equivalent to a change in the density. 

Let us first discuss the pair term, which depends only on the intcrparticlc distance r: 

vP^"{r) = [2Xd - AX^r^] V{r) + 2X\^ [V{r)]^ , (12) 

where we have used explicitly the (convenient) exponential form of the potential. In Fig. 1 wc show the form of this 
two-body potential. 

The form of the pair potential in the large A limit is simple: 

• For r > 1, the pair potential is zero because the exponential terms go very fast to zero. 

• For r < 1, the pair potential goes to infinity. Although the second term is negative, the third one is much larger 
than the first two and dominates. 

• Very close to r = 1, within a window that shrinks to zero in the large A limit, the potential becomes negative 
and goes to minus infinity. Indeed, for r = 1, its value is 2d\ — 2A^. The potential remains attractive (negative) 
on an interval of size 1/A. This attractive part of the potential is responsible for the fact that ground state 
wavefunction does not vanish (or is discontinuous) at the contact between hard spheres. Particles can stay 
close to contact paying a lot of kinetic energy, ©(A^), and gaining a lot of potential energy, also 0{X^). 

Let us now consider the 3-body term. It can be written as: 

^3-body(^^ x') = A2 X • x' V{r) V{r') , (13) 

where for a given triplet of particles i, j,j' one has x = x^ — Xj and x' — x^ — x^/, see Eq. (11). The form of the 3-body 
potential in the large A limit is simple: it is non-zero only if both particles j and j' have a non-zero overlap with 
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FIG. 1: Form of the pair potential v^'^"{r) in the mapped quantum problem that derives from the classical potential V{r) = 
Vo exp( — A[(r/cr)^ — 1]). (Notice that we have set Vb = 1 and cr = 1 in the text.) 

particle i. When there is a finite overlap (in the large A limit) the contribution coming from the 3-body term, which 
can be actually positive or negative, is always smaller than the one coming from the sum of the pair contributions ij 
and if . Thus, it can be neglected. The outcome coming from the pair interactions is that all particle configurations 
with finite particle overlaps are simply excluded by the Hilbert space. When both r and r' are very close to 1, of 
the order 1/A, the 3-body term will give rise to a non-negligible contribution similar to the one studied for the pair 
contribution. The main difference is that it can be attractive or repulsive depending on the relative orientation of x 
and x'. As a conclusion the interaction between particles is the hard sphere one plus a contact term that is sticky for 
the pair contribution and that can be repulsive or attractive for the triplet term depending on the geometry of the 
triplet. 

IV. SUPERFLUID, SUPERCRYSTAL AND SUPERCLASS PHASES 

Here we make use of the relationship between the classical and quantum models, and study the zero temperature 
phases of the quantum system introduced in the previous section III.D. We will work out the properties of the quantum 
phases by using the classical Gibbs measure defined by the square of the Jastrow wavefunction^^'^'*"'*^. 

We shall focus on a system of N identical bosons with the particular interaction, discussed in the previous section, 
that corresponds to a Jastrow state with an hard sphere potential. Having fixed a = 1, the control parameter in 
this case is the particle density p = N/V of the particles, or the packing fraction cj) — np/G. Let us recall first the 
phase diagram of classical hard spheres^^'^"'"'*'^ (see top of Fig. 2). At low density, the system is liquid, and upon 
increasing the density it undergoes a first order phase transition to a crystalline state, that is arranged in a face 
centered cubic (FCC) lattice. Moreover, a metastable glassy phase can be obtained in the classical problem if the 
density is increased fast enough or in presence of small bi-dispersity (binary mixtures). This glass phase can be 
compressed until the random close packing packing (RCP) fraction ~ 0.64. These three classical phases have their 
corresponding counterparts in the mapped quantum system. As we shall show, these are superfluid, supercrystal and 
superglass phases. The classical phase diagram and its quantum counterpart are shown in Fig. 2. 

Below we discuss the property of these three quantum phases in some detail with particular emphasis on their 
one-particle density matrix and Bose condensate fraction. We will treat the glass phase as a true equilibrium phase 
because, as we have already discussed, it will emerge when crystallization is avoided either by fast compression of the 
liquid state or considering suitable binary mixtures. 



A. Off-diagonal long range order and classical correlation functions 



Let us first recall general results on off-diagonal long range order (ODLRO)"**'^^. In particular, we shall show that 
for Jastrow-type wavefunctions there is a simple relation"''*'''^''** between the structure factor of a classical system 
and the one particle density matrix. This relation has a particularly simple form in the case of hard spheres. The 
condensate fraction is related to the amount of off-diagonal long range order in the quantum problem and can be 
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Glass transition 
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FIG. 2: Classical phase diagram for hard sphere systems (top) and its quantum counterpart (bottom). 



therefore related to an observable for the classical corresponding liquid. In appendix A 2 we generalize the following 
discussion of ODLRO to the case of binary mixtures. 

The one particle off-diagonal density matrix reads in term of the ground state wavcfunction "^q [V is the volume): 



7^(x,x') ^V I dx2...dxjv *g(x,X2, ...,X7v) *g(x',X2, ...,xjv) = V^'^-i V'j(x)^*(x') 



(14) 



where the rii and "02 (x) are, respectively, eigenvalues and (orthonormal) eigenvectors of the one particle density matrix 
7?.(x, x'), / dx' 7?.(x, x') V'i(x') = Tii ipi{'x.), see e.g. Ref 49. One can interpret rii as the fraction of particles in state i, 
and ^j Tii = 1. The largest eigenvalue no gives directly the condensate fraction. 

If ^G is given by a Jastrow wavefunction with a hard sphere potential, this expression is particularly simple in 
terms of a classical equilibrium correlation functions: 



7^(x,x') = 



V 



N{N+l) Z 



"^+1 y(x-x')/T 



N 



5](<5(x-x,)5(x'-x,)) = 



^Af+l 



ij^j 



ZNiN + 1) 



p-V(x)p(x')2/(x,x'), (15) 



where p(x) = (X]i^(^~ x;^)) and j/(x, x') = e^''^~^ ^/-^^(x, x'), and g(x, x') is the usual pair correlation function^". 
Notice that this expression follows by considering a system of A^-|- 1 particles, with two of them fixed at x and x', and 
that the difference between 2/(x, x') and (7(x, x') arises because 7?.(x, x') does not vanish when x and x' are less than 
a particle diameter away. For hard spheres, there is a simplification and .g(x, x') = 2/(x, x') for |x — x'| > 1; however, 
for |x — x'l < 1 the pair correlation function vanishes while j/(x, x') is finite^". 

The prefactor in Eq. (15) is related to the fugacity of the system and it is easily evaluated in the large N limit: 



^N 



TV! e^scw/y) 



^N+l 



Zn(N + 



^ e^ipSiP)] = ^-ie^[pSc.(p)] ^ p-if^p) ^ 1 



(P) 



(16) 



where S{p) = 1 — logp + S^xip) is the entropy of classical hard spheres as a function of the average density p — N/V, 
S'cx(/o) is the excess entropy with respect to the ideal gas, and z{p) is the fugacity. Using this relation we find the 
final, quite simple, expression: 



7^(x,x') = .f(p)^^^^^y(x,x'). 
P 

For a homogeneous phase, p(x) = p and 7?.(x, x') becomes a function of |x — x'| only, 

7^(|x-x'|)=/(p)y(|x-x'|) 



(17) 



(18) 
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FIG. 3: Condensate fraction as a function of the packing fraction for the hard sphere Jastrow wavefunction. Full (black) curve: 
superfluid phase, Eq. (20). Dot-dashed (blue) curve: supercrystal phase, Eq. (24). Dashed (red) curve: supergiass phase, 
Eq. (39); note that the theory slightly overestimates the random close packing density with respect to the commonly accepted 
value of = 0.64. The blue dot marks the location of the melting transition, where liquid-crystal phase coexistence begins, 
while the red dot marks the glass transition density. Inset: the enhancement factor no//(p) due to the spatial inhomogeneity 
of the crystal and glass phases. Note that this factor is quite large, of the order of lO'^ — 10^ at the melting and glass transition 
densities. Still, the condensate fractions we find are extremely small, probably due to the classical-like nature of solids with 
the Jastrow hard sphere wavefunction. 

and the condensate fraction is just lim|x_x'|^oo 7?-(|x — x'|) = /(p), since 2/(x, x') -^ 1 (like the pair correlation 
function) as |x — x'| -^ oo. This expression shows that as long as the fugacity is finite, i.e. the pressure is finite, we 
shall find a non-zero condensate fraction^*. As a consequence, we expect a superfluid phase at low densities, then a 
first order transition into a supersolid crystalline phase as the density is increased. In the case of non-identical hard 
spheres, e.g. bi-disperse, or for extremely rapid quenches, the system should end up in a supergiass phase at high 
density, and the condensate fraction will vanish only at close packing. In the following, we will apply the results above 
and study in detail the superfluid, supercrystal and supergiass phase. 



B. Superfluid phase 



At low packing fraction the classical system is in the liquid phase. Analogously, the quantum system is in a quantum 
liquid phase. As discussed above, we expect it to be superfluid. In the following we shall compute its corresponding 
condensate fraction. In order to obtain a heuristic semi-quantitative expression for the condensate fraction, we used 
the Carnahan-Starling approximation^'' for the classical hard sphere system. This is a kind of virial resummation 
known to work very well for describing the properties of the hard sphere system in the whole fluid range. Within this 
approximation one finds^'' 



pS^Ap) 



1 



log (Zjv/y 



N\ 



60 0(4 - 30) 



Then Eq. (16) leads to the following expression for the condensate fraction uq: 

(-8 + 90-302)" 



na = f{p) = exp 



(1 



This is plotted in Fig. 3. Note that the freezing crystallization transition takes place at a packing fraction 
for hard spheres, see Fig. 2. 



(19) 

(20) 
0.494 



C. Crystal phase 



For > 0.545, the classical hard sphere system is in the crystal phase (for 0.494 < < 0.545 there is phase 
coexistence between the liquid and the crystal). Correspondingly, the quantum system displays a quantum crystal 



phase. In the following we shall compute its condensate fraction. The crystal phase is inhomogeneous but still 
y(x, x') -^ 1 at large |x — x'|; therefore we have^^ 

hm 7^(x,x') - f{p) ^M^ ^ noVM^)ro{^') , (21) 



and from this factorization one obtains that y/ng V ipo{x.) = yj f{p) p{x.)/p (notice that the eigenvector V'o(x) is real 
valued, as 7?.(x, x') is real and symmetric). It thus follows from the normalization of "00 (^) that 



11/" 

"0 = f{p) —y dxp{x) 



(22) 

A rough estimate of no is given by f(p) since f(p) < uq < (maxy ^^^ ) /(p)- The first inequality follows from 

p^ = [y J dx p(x)l < T7 / dx p(x) . However, because of the inhomogeneity of the density profile, the contribution 
to no coming from the integral of the density squared may be quite large. In order to estimate it we assume that 



where Ri are the FCC lattice sites and the Gaussians are roughly non-overlapping. It has been shown by numerical 
simulations that this approximation is rather good^^ and values of A have been computed by density functional 
theory^^. Using these results one finds 

..^/W^^/oix„x)=^M_L_^^ (24) 

with f{p) given in (16). 

To compute /(p) we used the phenomenological equation of state for the FCC crystal phase of classical hard spheres 
proposed by Speedy^"^: 

P ^ z-b 

— = -. a , z = (p/cpFCC , (25) 

p 1 — z z — c 

where (fipcc — 7r/(3\/2) — 0.74 . . ., a — 0.5921, b — 0.7072, c — 0.601. The entropy can be obtained by integrating 
this relation with respect to density as detailed in Ref. 53 where the integration constant is also reported. Values of 
A have been taken from table I of Ref. 52. They have been fitted to a polynomial to get 

^TA = fli ((/)FCC - 0) + a2 (0FCC - 0)^ + as (0fcc - 0)^ , (26) 

with ai = 0.35, 02 = —0.866 and as = 3.58. The final result for no is reported in Fig. 3. 

It is worth to note that the crystal phase is not commensurate, except at close packing. Therefore our results are 
consistent with the general statement by Prokofev and Svistunov, that commensurate crystals are not supersolids^*. 

D. Glass phase 

In order to understand the superglass phase, one must first translate the results of the classical hard sphere problem 
into the quantum one. Thus, let us start by discussing the classical results. 

It is well known that the classical hard sphere problem may undergo a glass transition at a packing fraction (J)k ^ 0.6 
if compressed sufficiently fast^^'^^""^^''*^ or in case of mixture of different particles^^''^^. This has been found analytically 
under some approximations^^'"^"'^^"'^^, in simulations'^^'^'* and in experiments on colloidal systems^^'^^. The classical 
phase diagram is reported in Fig. 2. 

The physical mechanism behind this transition in finite dimension is still unclear. However, approximated theories 
able to provide good quantitative predictions have been developed. The most successful is the so-called Mode- 
Coupling Theory (MCT)'^"''*^'^^ that correctly describes the dynamics of the Langevin system Eq. (1) for densities 
slightly smaller than the glass transition density. A complementary approach is based on the replica method^^^^'*'^* 
and gives predictions on static observables at and above the glass transition density. We will discuss in the following 
what predictions can be derived from these theories for the quantum problem. 
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Our strategy will consist in deriving results for the quantum glassy phase starting from the known results for the 
classical glassy phase. We will also provide quantitative results using the replica approach which has been shown to 
give a reasonable description of numerical simulation results'^^'^'*. Our working hypothesis is that an ideal classical 
glass transition indeed takes place: at low density (or packing fraction) the system is in the liquid phase and above a 
critical density (assuming that crystallization is avoided) there is a thermodynamic transition toward an amorphous 
state of matter which is the glass state. In the glassy state the system can be frozen in many different amorphous 
configuration or states. Note that assuming that the classical glass transition is a real thermodynamic phase transition 
just simplifies the presentation. In fact, even if this is not the case in reality, one can translate and generalize all the 
following discussion: amorphous thermodynamic states will then become just metastable amorphous configurations 
in which the system is trapped on the relevant experimental timescales. 

1. The decomposition of the Gibbs measure in pure states 

Following the strategy outlined above, we will therefore assume that at the (classical) glass transition, the Gibbs 
measure Pg({x}) = |\I'g({x})P splits into a very large number of thermodynamic states Pq,({x}): 

|vl/G({x})p = ^^„„P4{x}). (27) 

a 

An operative description of Pq,({x}) is the Boltzmann measure obtained by coupling the system to an infinitesimal 
non- homogeneous external potential, Va, that forces the system into the state a: 

P„({x}) = hm , Z„ - hm / d{x}e-^"/^-^^°/^ . (28) 

Note that clearly the potential Va must be chosen symmetric under particle exchanges. From the above definition it 
follows that the states Pa are also completely Bose symmetric like the full Boltzmann measure Pq- 

The tricky aspect of glassy systems is that this external potential (or field) is not known a priori as it is for, say, the 
ferromagnetic case. The reason is that it is random as the state it pins. Wa is the thermodynamic weight of the state 
a. It equals Za/Zf^f, where Za is the partition function in the presence of an infinitesimal pinning field forcing the 
system into the state a. Zjq is instead, as before, the partition function in the absence of any pinning field. Hence, 
it contains the contribution of all different states: Zj^ = ^^ Za and then ^^ Wa = 1. 

Since the Pa({x})s are equilibrium steady state distributions one can obtain an eigenfunction with minimum (zero) 
energy using the similarity transformation introduced in Eq. (9): 

v|/„({x}) = Na y^ e^"/2T p^(|^^) ^ ^^ P„({x})/v1/g({x}) . (29) 

J^a is a normalization constant that has to be fixed by imposing that ^q({x}) is normalized: 

/„~2UN/T-2eVa/T 7 r 7 

d{x} Z^ e^"/^ '- ^ Ml Inn ^ j d{x} e-^«/^--^"/^ = A^^^ 

Note that we have used explicitly that one always recovers Za whether one uses e or 2e. The outcome is that 
Ma = ^J Za/Zjsi = ^Wa- Usiug this rcsult and plugging the expression (29) into (27) one finds: 

^G{{^}f = *g({x}) Y. V^ Vl/^lx}) . 

OL 

By factoring out ^g({^}) one finally finds^*^: 

v^G^xD-J^V^vl/^dx}). (30) 

The ground state wavefunction is therefore the coherent sum of the wavefunctions corresponding to states a. 

However, no interference is present in the thermodynamic limit and all cross-terms in the square of the wavefunction 
can be dropped. This can be shown as follows: first, by plugging the expression of Ma into Eq. (29) one finds that 
^a({x}) = ^Pq,({x}). Using this result and noticing that the square of the previous equation (30) has to give back 
Eq. (27) we find that the interference (cross-products) terms have to be zero. This result can be understood in a 
simple heuristic way: configurations on which the wavefunction ^q, is concentrated have extremely small weight in 
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any other ^/j. Any given interference term a, /3 is expected to give a contribution decreasing exponentially fast in 
iV. Since the number of states, hence the number of couples a, (3, increases at most exponentially in N'^'^~^^/'^ (see 
Rcf. 59) one finds that the interference terms can be neglected in the thermodynamic limit. 

Thus, the probability of finding the quantum system in a state a is not different from that of the classical problem: 
both are given by Wa in absence of a pinning external potential. The pinning potential Va introduced in Eq. (28) 
induces a corresponding pinning potential Va via the mapping introduced in Eq. (8). The effect of the pinning 
potential is to concentrate the Jastrow wavefunction on state a; equivalently, in the quantum problem, the pinning 
potential Va lifts the degeneracy between quantum ground states and selects ^^ as the unique ground state of the 
system. As a consequence each one of the ^^s has to be interpreted as a possible state of the system (a more direct 
dynamical interpretation that shows that the system does not escape by tunnelling once it is in a state a is presented 
in Sec. V) . The glass phase in the quantum system is thus a random solid in the same sense as its classical counterpart. 
However, much as in the case of the supercrystal, ODLRO can still develop, leading to a non-vanishing condensate 
fraction, which we study below. 

Translating the classical results for hard sphercs"^*^ in the quantum case we thus find that the quantum ground 
state is unique for (j) < <J)k and is degenerate for (j) > (J)k but the logarithm of the number of ground states is 
subextensive and therefore the entropy remains zero above 4>k- It follows that the phase transition is not manifested 
by a non-analyticity in the free energy of the system. However, suitable correlation functions (such as the "point to 
set" or "dynamical" correlation functions^°'^^) should display a growing correlation length at the transition. Also 
the structure function g{r) shows a (weak) discontinuity at (p^^ but it does not display any long-range order. In the 
following we will study the properties of this glassy quantum state, in particular the ones related to ODLRO. 

2. Condensate fraction and corresponding Edwards-Anderson parameter 

Using the results above on the ground-state wavefunction one can readily obtain the one particle off-diagonal density 
matrix. As interference terms can be neglected, it reads: 

TZ{x,X.') ^^Wa 7^a(x,x') == V^Wa^Ui^a 1pi,a{x) i'lai^') , (31) 

a a i 

where the rii^a and the ipi^ai'^) are, respectively, eigenvalues and eigenvectors of the one particle density matrix 
7?,q(x, x'). (Because the wavefunctions ^q, are real, the 7?.q,(x, x') are real and symmetric, and the eigenvalues and 
eigenvectors are real.) The n^.Q have a physical interpretation: they are the fraction of particles in the quantum state 
i within the thermodynamic state a, and ^^ ^ Wa ni_a — 1 accounts for all the particles in all possible states a with 
probabilities Wa- 

For large separations |x — x'| ^ oo, one has for each TZa 

, hm 7^„(x,x') = ./a(p)^^^^^%^-r^o,„FV'o,a(x)V'oa(x'). (32) 

|x— x'|— >oo p ' 



Proceeding similarly to the case of the crystal, we obtain i/'o,a(x) = \iv Pa{^)/P- We use the normalization of 
the eigenvector -00,0: to write 

1 1 f 2 

no,a = fa{p) —y dxpa{x) , (33) 

and for the total condensate fraction we get 

1 1 r 

no = —y^ Wa fa{p) 77 / dx pa{xf . (34) 

Pa ^ J 

We expect that the fa{p)s are independent of the state a, since they are thermodynamic quantities related to the 
addition of one single particle to the system [Recall Eq. (16)]. In order to quantify the inhomogeneity of the glass 
state one can introduce the correlation function of the density and condensate wavefunction fluctuations that read 
respectively: 

g{x,x') = p-'^^WaPa{y^)Pa{y^') , G^ (x, x') = ^ Wa-^O.a (x)-0O,a (x') • (35) 



12 



100 



10 r 



g(r): 



0.1 



r 


' 




- 


: 


y 


— 9=0.63 

— 9=0.6635 

— 9=0.6905 


- 


1 







1000 




FIG. 4: Correlation function of the density (left) and of the condensate wavefunction (right) for the quantum system corre- 
sponding to classical hard spheres at different densities in the glass phase. Details of the computations are in Ref. 33. Note 
that these quantities can be measured in Quantum Monte Carlo numerical simulations, see e.g. Ref. 20. 



Actually, using eq. (32) we find that they are related one to the other by a simple proportionality relation: 

5(x,x') = -— — V'fi;a?/;o,a(x)?/;o,Q(x') = -— — G'^(x, x') . (36) 



Substituting these definitions in Eq. (34) we get 

= /(P)y / rfx— ^wa Pai-x-f = f{p)y j dx g(x,x) = f{p)g{0) , 



no 



(37) 



where we used that g(x, x') = 5(|x — x'|) — g{r), as translational invariance is restored after averaging over all possible 
states. Note that g(x,x') plays the role of an Edwards- Anderson order parameter for the glass state. For instance 
one can define 



9ea = X! ^" 17 / ^^ 



Pa(x) 
P 



-1 =5(0)-l 



(38) 



Since the density and condensate wavefunction fiuctuations are proportional g(x, x') also represents (up to a propor- 
tionality constant) the inhomogeneity of the condensate wavefunction. It is a quantitative measure of how much the 
condensate is amorphous. It can be computed in numerical simulation. Look for example to Fig. 3 of Ref. 20 which 
provides a visual representation of the inhomogeneity captured by ^(x, x'). 

In order to obtain quantitative results on the condensate fraction we have to compute g{0) and f{p). The replica 
method'^^, in particular the small cage expansion^*, has been successfully applied to describe the glassy phase of hard 
spheres'^'^'^''. These two quantities can be indeed obtained using the procedure detailed in Refs. 33,34. Since this is 
well documented^'* we do not reproduce the computation and just quote the final result which is very similar to the 
one for the crystal: 



no 



fip) 



p (47rA)3/2 ' 



(39) 



where A is the so called cage radius and is a measure of particle vibrations. Hence, by using the equation of state for 
the glass and the values of A reported in Ref. 33,34 we finally obtain a quantitative result for the glass condensate 
fraction which is reported in Fig. 3. Note that the results of Ref. 33,34 depend slightly on the approximation that is 
used to describe the liquid; for consistency, we used the results obtained using the Carnahan-Starling approximation 
as we did for the liquid. 

The replica method allows to compute both g{r) and 'g{r) in the glass phase, but within a different approximation 
scheme known as HyperNettcd Chain (HNC) approximation^". The results, taken from Refs.^^'^'^, are reproduced in 
Fig. 4. 
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V. QUANTUM SLOW DYNAMICS AND THE APPROACH TO THE QUANTUM GLASS 

TRANSITION 

In the following we shall focus on the real time dynamics of the superfluid phase when the transition to the superglass 
is approached. Also in this case, the knowledge of dynamical correlation functions for Brownian hard spheres will 
allow us to obtain results on dynamical correlation functions of the corresponding quantum problem. In the next 
section we will explain how the mapping works for dynamic observables. 

A. Mapping from Brownian dynamics of hard spheres to real time dynamics of the quantum model 

In order to show how one can obtain information about the real-time quantum dynamics from the Langevin dynamics 
(1) it is useful to introduce a bracket notation for the Fokker-Planck problem; we define P{{x},t) = ({x}|P(i)), and 
we denote by \G) the Gibbs distribution (5), Pcdx},*) = ({x}|G) = exp[-C/A,({x})/T]/Zjv, such that Hfp\G) = 0. 
We also denote by (+| the constant state, (+|{x}) = 1. Note that 

hIp ^ e^"/^ Hfp e-^"/^ , (40) 

which is consistent with Eq. (6) and the fact that H is Hermitian. Then 

Hfp\G)=0 => 0={G\Hlp = {G\e^"/^ Hppe-^"^^ ^Z]^\+\HFpe^^"^^ ^ = {+\Hfp (41) 

i.e. (+1 is a left eigenvector of Hpp with zero eigenvalue. 

One observable that is particularly interesting to characterize the quantum dynamics is the dynamical structure 
factor, FQ(q,t) which is the time dependent correlation function of a Fourier component of the density operator, 
/9q({x}) — ^i e^"^'^'. We shall show that the dynamical structure factor for the quantum problem is related to the 
imaginary time analytic continuation of the dynamical structure factor, Fci{q,t), for Brownian hard spheres. 

Observing that ({x}|e^*^^^ |{y}) is the probability of going from {y} to {x} in time t, we can write the Brownian 
correlation function as follows: 

r „ p-/3'7„({y}) 

F,,(g,t) - {P,{t) P-M) ^ J d{^}d{y} p,({x}) ({x}|e-*^--|{y}) p_,({y}) — 

= (+1 P, e-'"-- p^,\G) = (+|e*«- p, e-*^- p_,|G) = (+|e-^"/(2^) e*« p, e''" e^^^^'^^ p_,|G) (42) 
= (0|e*^ p, e-'" p_,|0) = J2mp,\n)\' e-'i^.-Eo) = / ^ p,(c.) e^- , 

n ^ 

where |0) = y/Zjy e'^^/'^^^-'jC) is the quantum ground state (10), \n) are the excited states, and 

p,{uj) = 2^^ |(0|pg|n)|2 5{E„ -Eo-Lu) (43) 

n 

is the distribution of classical (inverse) relaxation times for the density fluctuations. 



We are interested in the quantum correlations; the quantum response function is given by' 



62 



n)P 



^-it(E„~Eo) _ ^it{E„-Eo) 



(44) 



Rgiq, t) = i9{t) (0| [e"« p, e-^*«, p-,] |0) = i9{t) ^ |(0|p,| 

n 

= -29{t) / -^ pq{uj) sin(wt) . 
It follows that the imaginary part of its Fourier transform reads: 

i?^(g,a>)-isgn(^)p,(|a.|), (45) 

and, using the quantum (bosonic) fluctuation-dissipation theorem at zero temperature, we get the quantum correlation 
function: 

SQ{q,uj) = \sgniu;)R'^{q,co) ^ ]- p,i\Lo\) ^ FQ{q,t) = T ^ SQ{q,uj) e^ * = H ^ P,{oj) cos(^t) (46) 

2^2 J_^ 2TT Jq 2TT 

This final simple expression allows us to obtain results on the dynamical structure factor starting from the distri- 
bution of (inverse) relaxation times. 
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FIG. 5: (Left) Long-time shape of the correlation function Eq. (49) Fci{q,t) ~ exp(— ^/t/ro.) and FQ{q,t) ~ 
exp( — -\/i/2Ta,) cos{^/t/2Ta) as a function of t/ra in a linear-log scale. (Right) Imaginary part of the response function of 
the quantum problem, Eqs. (48) and (49), as a function of frequency for different values of Tq (in arbitrary units). 



B. Results for the time dependent density-density correlator 

Assume that the density correlator has been normahzed in such a way that F{q,t = 0) = 1; then f^ ^pq{ijj) — 1. 
In the glass transition literature, different phenonienological expressions have been used in order to describe the 
slowing down of the relaxation of F{q, t) on approaching the glass transition. 

In dielectric spectroscopy one often writes 



/oo 
dlnr G,(lnT)e-*/^ G,(lnT) = 

OO 



2ttt 



(47) 



and to describe the slow relaxation (low- frequency) part it is assumed that Gg(lnr) contains a single time scale Ta{q), 
GqihiT) = {T/Ta{q)) Q (j / Ta(q)) ■ The scaling function g{x) might depend weakly on q. Then 



Fciiq.t) 



Fqiq^t) 



dxg(x)e-*/(-^"(?))=/,Ki/^a(g)) 



Kil^co)^- ^ 



dluT Gg(lnT) cos(i/r) 
1 



dx gix) cos . 

\XTa{q) 



fQ{t/ra{q)) 



9 



LU UJTa{q) \UJTa{q) , 

Several forms have been used successfully in the literature, see Ref. 63 for a review: 
1. Stretched exponential (3 = 1/2 - The simplest non-trivial case is fciiy) — e^^; in this case 



gix) 



/Attx 



fQ{y)^e-^y/^ COS \ 



V2] 



^g^-s-(xs)^cos{-.i3) sin[(xs)'^ sin(7r/3)] , 



2. Stretched exponential - For fciiy) — e ^ and generic /3 we have 

9{x) = — ds e-^-("^'' 

and in this case fgijj) must be computed numerically. 

3. Cole-Davidson - Another common expression, corresponding to 

sin(7r7) f x ^'^' 
9[x) - 



<x < 1 



IT \1 — X ^ 

In this case fgiy) can be computed in terms of Hypergeometric functions. 



(48) 



(49) 



(50) 



(51) 
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Here we limit ourselves to illustrate the qualitative behavior of Fq {q, t) on approaching the glass transition. For 
this, we consider solely the first case. (We checked numerically that the other forms give qualitatively similar results 
- oscillations in Fq are much more pronounced if one uses the Cole-Davidson form due to the frequency cutoff.) Some 
examples constructed using Eq. (49) are reported in figure 5. The classical correlation function decays from 1 to 
over the time scale Tq; the corresponding (real-time) quantum correlation function, directly accessible in numerical 
simulations, displays the same slow decay modulated by oscillations according to (49). The imaginary part of the 
response function, directly accessible in experiments, shows a peak at w ~ I/tq, whose amplitude increases oc Tq, upon 
increasing Tq. Note that usually in classical glassy systems these low- frequency features are accompanied by faster 
relaxations, related to "intra-cage" motion and/or fast molecular relaxations. These have been neglected here and 
in figure 5 but we expect them to be present in the quantum case as well. They would appear as secondary peaks 
at higher frequency, and typically their density dependence is weak. The complete time-dependence of the density 
density correlation function would therefore be a first rapid relaxation to a plateau value and then a second one which 
is the one studied in detail in this section. At the glass transition this second relaxation does not take place anymore 
since the relaxation time diverges (or it is larger than any experimental timescale). As a consequence part of the 
density fluctuations becomes frozen. The plateau value in the correlation function measure precisely that. Actually, 
the plateau value for a wave- vector q will coincide with the Fourier transform of the function §^(x — x') defined in 
Eq. (35). 

C. Condensate Fluctuations 

To conclude this section we will focus on the dynamical fluctuations of the condensate wavefunction. As we already 
discussed, in the glass phase the system can be in many different states, each one characterized by a different density 
profile and by a condensate wavefunction, '0o,q(x), that in our model is simply proportional to the classical density 
profile of the corresponding state. In Eq. (36) we defined the correlation function of the fluctuations of the condensate 
wavefunction when the system is frozen in a given amorphous state. This has the form reported in figure 4 and can 
be accessed by a direct computation of the condensate wavefunction as done in Ref. 20. It is important to notice that 
it can also be obtained from purely dynamical measurements as we will now show. 

Consider the time-dependent one particle density matrix; in second quantization it is defined as 7?.Q(t;x, x') = 
(V'^(x,i) ^(x',0)), where ^(x,t) is the standard bosonic annihilation operator. Within first quantization, it can be 
written as follows: 

7?.Q(i;x,x') ^V dx2- ••dxAT dx.2 ■ ■ ■ dx'^ *g(x,X2,- • • ,xn) (x2,- •• ,XAr|e"**-^|x2, • •• ,x'^) *g(x',X2, • • • ,x^) . 

(52) 
Clearly for t = Eq. (52) gives back Eq. (14). In order to understand the time dependent behavior it is useful to 
consider the evolution in imaginary time and use again the mapping on the Langcvin dynamics. From (6) we have 

(X2,--- ,xjv|e-*«|x'2,--- ,x'^) = [*g(x2,--- ,XAr)]-i(x2,--- ,XAr|e-*^--|x;,,... ,X^) *G (x'2 , • • • , x'^) , (53) 

where '^0(^2, • • • , ^n) is the Jastrow state (10) for A^ — 1 particles. Plugging (53) in (52) and using the explicit form 
(10), we finally obtain 



7^(t;x,x') = y /"dx2---dxArdx^---rfx'^ e-3 2:"2V'(x-x,) (x2, • • • ,XAr|e~ 



|X2i y^N/ 



X e"2 .^ 



.-hT.^=^v(^'-^[) e 



(54) 



■^N 



where Un-i is the interaction potential of the A^ — 1 particles and we have used the notation TZ for the imaginary 
time continuation of TZq . 

In the special case of a Jastrow hard sphere wavefunction, the factors of ^ are irrelevant. Then Eq. (54) becomes 

/^ p — C^Af (x',X2,-" jX^) 
dx2---dxArdx^...dx'^e-2:,=2^(x-x,) ^X2,... , xjy le"*^^^ |x^ , • • • ,x^) , (55) 

and has a straightforward interpretation in terms of Brownian dynamics of hard spheres: one should pick up a 
configuration {x'} from the equilibrium distribution, such that particle 1 is in x.'^ = x'. Then particle 1 must be 
removed, and particles 2, • • • , A" evolved according to the Langevin dynamics in absence of particle 1. Finally, one 
should attempt to reintroduce particle 1 at position xi = x. The function 7^(t;x, x') is the probability that the 
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attempt is successful, or in other words that there is a void at time t around x large enough to allow the reinsertion 
of particle 1. 

Close to the glass transition, as we already discussed, there is a huge separation of time scales in the classical 
dynamics between a "fast" intra-state relaxation and a "slow" relaxation corresponding to hopping between different 
states and characterized by a growing time scale r^. If we remove particle 1 at i = 0, the perturbation of the density 
field will not relax inside the initial state until i ^ Tq, where the state will change. This implies that 7^(i;x, x') will 
have a plateau at times Tfast <C i ;$ Tq, corresponding to the stationary probability inside the initial state at t ^ 0. 
This is given by 

^^(Tfast < t < Ta;x,x') CX ^ Wq V'O^qCx) V'O^qIx') OC ^(x - x') , (56) 

a 

where one has to average over all possible initial states. As obtained in the previous section for the density-density 
correlation function, we expect that in presence of a huge separation of timescales, as it is the case close to the 
glass transition (or in the glass state where the second relaxation does not take place anymore), one finds that 
'R-irtust <C i ^ Tq; X, x') ^ T^q (Tfast '^ t < Ta;:x-, x'). One way to understand this result consists in expressing TZ{t) in 
terms of the Fourier transform of TZq: TZ{t) = / 57 T^qi^) exp(— wt). If TZ displays a very long plateau this means 
that TZq{uj) contains two distinct contributions corresponding to w oc 1/rfast and lo cx 1/tq,. As a consequence, TZq on 
times intermediate between fast and slow timescales will coincide with TZ since the contribution from large frequencies 
(of the order or l/rfast) will have died out and the contribution from the very long frequencies (of the order of 1/tq,) 
will be the same since exp(— iwt) ~ exp(— wt) ~ 1. As a conclusion, TZq will display a plateau whose extension will 
become infinite beyond the glass transition. The value of TZq on the plateau corresponds to the fraction of frozen 
condensate wavefunction fluctuations and equals G^ (and hence is proportional to g). Therefore this quantity can 
also be computed in a dynamic framework without introducing replicas. 

VI. SUPERFLUID PROPERTIES AND PHASE DIAGRAM OF REALISTIC SUPERCLASS PHASES 

The conclusion from the previous sections is that the ground state of the model can be a liquid, a crystal or a 
glass and all these phases are characterized by a finite condensate fraction no. However, the study of the superfiuid 
properties requires also the knowledge of excited states, or at least of the excitation spectrum. The latter is related 
to superfiuid properties by the celebrated Landau argument^'* that predicts for the critical velocity 

Vc < min[e(fc)/fc] . (57) 

k 

In the zero-temperature liquid He** phase, the excitation spectrum is linear at small k and has a minimum at larger 
k, therefore Vc is finite and the system is superfiuid. 

Unfortunately, for the Jastrow wavefunction discussed above, Eq.(lO), one can show that the excitation spectrum 
is quadratic at small fc, at least if the potential V^(x) has finite integral. Therefore Eq.(57) implies that Vc — Q and 
the system is not superfiuid, much as it happens for an ideal Bose gas despite a condensate fraction equal to 1 at zero 
temperature. 

The quadratic spectrum of Jastrow wavefunctions can be related to the following properties of the Hamiltonian (7): 

1. The ground state energy per particle e{p) of the Jastrow ground state is always zero, therefore the pressure 
P — p^^ = Q and the compressibility is infinite, Xr^ = Px" ^ ^' 

2. Consequently the sound velocity c = '^1 ^/pxr — 0, i.e. there are no phonons and the linear part of the spectrum 
at small k is absent; this is because no restoring force for density fluctuations is present if e(p) is independent 
of p. 

3. The static structure factor of the Jastrow wavefunction has the property linifc^o S{k) 7^ 0; therefore, the 
Feynman relation®* e(fc) = fc^/(2m5(fc)) gives e(fc) cx fc^, consistent with c — 0. 

4. Finally, it is possible to identify a hidden symmetry in the problem, related to the special form of the potential 
in (7), that is responsible for non-trivial cancellations in the Bogoliubov low-density perturbation theory for (7) 
around the ideal gas limit. Again these cancellations are responsible for the absence of the linear part of the 
spectrum. 

When using Jastrow wavefunctions as variational functions for liquid Helium, a classical strategy^^ to reintroduce 
phonons is to add a non-integrable part to the Jastrow potential, such that y(x) ^ |x|~^ at large |x|. In this way 
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FIG. 6: Schematic behavior of e(p), -P(pi), 
that correspond to the same density pK- 



fi{p) and P{p) across the glass transition. Note that there are values of P and p 



S{k) ~ fc at small k and the excitation spectrum is linear. This is also quite natural since long range correlations 
are expected in the ground state of generic quantum systems even if the original interaction is short ranged. The 
inclusion of these terms in our formalism is possible and their presence does not influence much the results for the 
properties of the glass, e.g. for g{r) and g{r) discussed in section IV D, except at large r/small k. 

However, the classical-quantum mapping will give in this case a quantum Hamiltonian (7) with very long ranged 
interaction, while we would like to keep the original local nature of the quantum Hamiltonian. Therefore we consider 
an alternative way to solve the problem: we propose to introduce a perturbation of the quantum Hamiltonian (7) by 
adding a small potential term AVAr({x}) = X]i<i Su{yii — Xj) and treat it in perturbation theory. 

The situation here is quite similar to Bogoliubov low-density perturbation theory^^: indeed, the ideal Bose gas has 
a finite condensate fraction (actually equal to 1) but is not superfluid because its excitation spectrum is quadratic. 
Once an infinitesimal interaction is added (or in the very dilute regime), the spectrum immediately becomes linear 
and the system becomes superfluid®^. Hence by analogy we argue that perturbation theory can be applied in our case. 
As a check of this argument we verified that in the low-density limit, as discussed above, the Bogoliubov theory of the 
unperturbed model leads to a quadratic spectrum due to a hidden symmetry. On the contrary, the perturbation AVn 
breaks the hidden symmetry, and in its presence the usual Bogoliubov theory applies and leads to a linear spectrum. 

Using perturbation theory at first order, it is straightforward to show that 



(58) 

One can show by explicit computation for a 
the ground state energy is continuous but its 



g(r) is the correlation function of the hard sphere liquid at density p. 
suitable specific form of 6u{r) that, at the glass transition density pK, 
first derivative has a jump. 

The pressure P{p) (as a function of the chemical potential p) is the Legendre transform of e{p): 

P{p) = max[/?/i - pe(p)] , 
p 

therefore 



(59) 






(60) 



From the above expressions one can see that -P(/x) is continuous and convex as e(/9), while p{p) and P{p) have a jump 
in pk, see Fig. 6. Therefore the glass transition, that is a second order transition in the classical case, looks like a first 
order transition in the zero-temperature quantum problem. This is similar to previous results on quantum mean field 
glass models with quenched disorder, see Refs. 66,67. On the other hand the properties of the first order transition 
are quite different. In the case of mean field quantum glass models the glass and liquid are really different phases and 
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the glass phase does not appear (via density fluctuations that are frozen on timescales diverging at the transition) in 
a continuous way from the liquid, contrary to what happens for the superglass. 
The sound velocity is now determined by 



2 _ 1 _dP _ d 
PXT dp dp 



c 



2 '^^ 

dp 



^0. (61) 



General arguments^^ show that if c 7^ 0, then S{k) < fc/(2?nc); assuming equality the Feynman formula gives e(k) — ck, 
consistently with the existence of sound waves. 

Therefore the system in this case has a finite critical velocity and is superfluid. Note that in the glass the first peak 
of the structure factor, that determines minfc[e(fc)/fc], is close to the one of the liquid, suggesting that the critical 
velocity should stay close to the one of the liquid. Still it is known (e.g. in Helium) that the precise determination 
of the critical velocity is complicated and depends on the geometry, therefore estimating precisely its value in the 
superglass phase is beyond the scope of this paper. 

VII. CONCLUSION 

In this paper we presented and analyzed a concrete model of interacting bosons that displays unambiguously, in 
addition to superfluid and supercrystal phases, metastable superfiuid and superglass phases. This shows concretely 
that a system can be at the same time glassy, i.e. displaying very slow dynamics for the structural degrees of freedom, 
and supersolid, i.e. showing a superfluid component. Note that there is no paradox, exactly as there is none for the 
supercrystal. Indeed, Leggett in his original paper'^ "Can a solid be superfiuid?" explicitly mentioned the possibility 
of super-amorphous solids. 

The exact ground state of the model we considered is known by construction: it is a Jastrow wavefunction with 
a hard sphere potential form. We constructed the quantum Hamiltonian for which this wavefunction is exact via a 
mapping to Brownian motion of classical hard spheres, following the connection between the Fokker-Planck operator 
for the classical stochastic dynamics and the Schrodingcr operator describing the associated quantum system. This 
mapping allows us to understand the physics of the superglass phase by using results from the well-studied problem 
of densely-packed hard spheres. 

Our findings, which were summarized in Sec. II, are most simply conveyed through Fig. 2. By changing the density 
of particles, one goes from a superfiuid phase to a supercrystal via a first order phase transition. If the density is 
increased fast or in case of (binary) mixtures, one reaches a metastable superfiuid state, and if the density is further 
increased, one reaches the superglass phase. 

The analysis we carried out for the hard sphere Jastrow wavefunction can be extended for other systems for which 
the Jastrow wavefunction corresponds to a classical potential leading to glassy dynamics. Classically, mixtures may 
become glasses at certain compositions, for example 80%/20% mixtures of two species of particles interacting via 
Lennard- Jones potentials^". The problem of binary mixtures is interesting in that one could potentially realize the 
system in cold atomic gases, where the relative fractions of two species could be controlled. We have worked out the 
correspondence between the classical and quantum problems for binary mixtures in Appendix A. Following the same 
reasoning as in our work on the hard sphere wavefunctions, the binary mixtures can display superglass behavior as 
long as the mapped classical problem does so. 

The usefulness of the classical to quantum mapping is that it allows us to make precise and concrete statements 
about the nature of the ground state of the devised Hamiltonian. By construction, the ground state energy is exactly 
zero, and thus one must add a small perturbation so as to obtain a non-zero speed of sound, as we have done 
in Sec. VI. We have added a 2-body potential as perturbation, because we could make use of the knowledge of 
the density-density correlation g{r) in obtaining the effect of the additional interaction on the ground state energy, 
pressure, compressibility, and sound velocity of the system. It is an interesting possibility that if, as opposed to 
the 2-body perturbation, one adds a 3-body potential, one could compensate the 3-body term that arises from the 
classical to quantum mapping. In other words, we speculate that by removing (or reducing) the 3-body term by a 
compensating perturbation, one could achieve a quantum Hamiltonian with a glassy ground state, only 2-body terms, 
and a finite sound velocity. This could perhaps be tested numerically, by considering the system with the potential 
in Eq. (12) alone. Moreover, if the short-ranged "sticky" part of the Eq. (12) potential is removed, this may not 
substantially change the results. At least, we expect the system to be glassy as the original one (the question of 
superfiuidity is instead a more tricky one and needs further investigations). The intuitive reason is as follows: for 
dense packings and a hard sphere potential, if one turns down the kinetic energy (say by "increasing" the mass of the 
particles), thus making the system more classical, one should recover the physics of the classical hard sphere system. 
Because the packings are dense, turning up the kinetic energy (say by "decreasing" the mass of the particles), should 
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not favor very much transitions between a glassy state a to another (3. Because of the hard sphere potential (basically 
infinite on the scales of the kinetic energy), we do not expect, at a dense packing, that quantum fluctuations will 
reduce considerably (at least for not too small m) the timescales to escape from classical glassy configurations. 

Concerning the superglass transition, we found it to be first order but quite unusual, ft is different from the quantum 
first order transitions found for mean field glassy models in Rcfs. 66,67. ft would be interesting to develop and study 
mean field models (or mean field approximations) able to reproduce our results and the superglass phase transitions. 
One of the main motivations is that they could allow one to analyze the finite temperature regime which is clearly 
out of reach of the approach developed in this work. In particular, it would be very interesting to apply Quantum 
Mode Coupling Theory^^ (QMCT) to the model we focus on. It would be very interesting to know whether QMCT 
predicts quantum correlation functions that are connected to the one predicted by classical MCT via the classical 
quantum-mapping we employed. If yes, this would mean that QMCT is able to capture this unusual first order phase 
transition toward the superglass phase and it would allow one to study its finite temperature extension. 

Finally, let us discuss preliminary results for He** which is one of the original motivations of our work. Jastrow 
wavefunctions have been used extensively to describe the ground state of He** in both fluid^^'^''' and solid''® phases. 
Therefore it is tempting to try to study the superglass phase of He^ by using the Jastrow wavefunction as a variational 
wavefunction. For He** it was found that the liquid phase is well described by a potential V{r) ^ r~^. We then 
computed the glass transition density for this potential following Rcf. 69. We found that the superglass transition 

takes place at a density pK ~ 0.45A and a very small jump in the first derivative of e(p) at pK, the relative variation 
of e'{p) being of the order of 10~^. Unfortunately, the value of pK is 10 times higher than the density at which the 
superglass has been observed in Rcf. 20. However, it is well known that the Jastrow wavefunction overestimates 
the liquid-crystal transition for He* and, more generally, provides a poor description of the solid phases. A more 
refined investigation should involve Shadow- like wavefunctions^^, and we plan to report on such a study in a future 
publication where more details on the Jastrow wavefunction will also be given. 
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APPENDIX A: QUANTUM-CLASSICAL MAPPING FOR A BINARY MIXTURE OF BOSONS 



Here we work out the details of the connection between the quantum and classical fluids for a mixture of two types 
of particles, A and B. We start again with the Jastrow type wavefunction written in terms of the interaction among 
N — Na + Nb particles, through a potential 

t/({x}; {X}) = ;7(xi, . . . ,XAr,;Xi, . . . , X^ J - ^ y^^(|x,-x,|) + ^ Vbb(|X,-X,;|) + ^ T/^b(|x,-Xj|) , (Al) 



i>j 



I>J 



where i,j = 1, . . . , Na and I, J — 1, . . . , Nb- The dynamics of a classical mixture of particles interacting through 
this potential is clearly a particular case of the general Langevin equation specified by Eqs. (1) and (2), where we set 
7i = J A and 77 = 7b , and the form of the Fokker-Planck operator follows straightforwardly from Eq. (4) . 
The quantum Hamiltonian is, from (7), 



^-E 



Pi 



2mA 



P\ 



E^+^(w={x})' 



2to 



(A2) 



where mA,B — f^lA.Bli'^T) and, setting T = 1 and h~l, 



V({x};{X}) = V(xi,...,XA,^;Xi, 



,X 



Ne 



^^■Z{\i^.uf-viu],J-^[\iS.vf-v',v] 



The ground state state of this Hamiltonian is 

vI/g({x};{X}) = 



\/ZnZn^ 



exp 



--t/({x};{X}) 



(A3) 
(A4) 
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FIG. 7: Form of the pair potential v^'^^^(r) in the mapped quantum problem that derives from the classical a = 6 and /? = 12 
Lennard- Jones potential, (e = 1 is set for simplicity.) The potential is repulsive at large and small distances, but to see the 
repulsion for large r one needs to zoom closer, as shown in the inset. 



The quantum potential V({x}; {X}) will again have contributions in the form of 2-body and 3-body terms: 



v = E 



AA pair 



Y^ vff P""' + J2 ^u"" "'" + 3-body terms . 



i>j I>J 

The pair potential is a simple generalization of Eq. (11) 

AA pair _ 



i,I 



V. 



y 



2^ {-^ ^AAin,) - VUini) + \ [v:,A{n,)? 



V 



BB pair 

IJ ^ 



1 



d-1 



1 



2m_B I rij 2 



AS pair mA+niB ( d-l , „ 1 , 2 

AniAmB I r,j/ 2 



''i/ 



(A5) 

(A6a) 
(A6b) 
(A6c) 



and similarly for the 3-body terms. 



Quantum model associated with Lennard-Jones binary mixtures 



Let us construct the quantum model associated to a Lennard-Jones binary mixture, for which the underlying 
classical system can be glassy, for instance certain 80%/20% mixtures of A/B particles™. Let us consider classical 
potentials of the form Vp^p^{r) = —ep^p^ [{r/crp^p^)~°' — {r/apj^P2)~^], where Pi,P2 = A or B, with three energy 
scales, eAA,£BB and eab, and three characteristic lengths, aAA^cTBB and ctabj in the problem. The exponents 
(3 > a > make the classical potential attractive at long distances and repulsive at short ones. The corresponding 
quantum pair potential is 



V 



Pi P-i pair 



(r) 



mp^ +mp2 ep^P^ 



,-2 

+£PlP2 



^ {{a{a + 2 - d) (r/ap,pJ-("+2) - (/3(/3 + 2 - d) (r/ap,pJ-(''+2) 



Y (r/ap,pJ-(^"+^) - aP (r/ap,pJ-('^+/'+^) + ^ (r/ap,pJ"(2''+2) 



(A7) 



Notice that the pair potential is repulsive for both small and large distances: t;^i-P2 pan-^^-j ^ f^^i^^^^^\^-{a+2) ^^^ 
rJGp^P^ > 1 and v^i ^2 pair (-j,-j ^ j-j,^^^^^_^-j- (2/3+2) j-q^. r/ap^p^ <^ 1. At intermediate distances, r/ap^p^ ~ 1, the 
potential can become attractive, as illustrated in Fig. 7 starting from an a = 6 and /3 = 12 classical Lennard-Jones 
potential. 



21 

2. Off-diagonal long range order in binary mixtures 

We can extend the Penrose-Onsager definition of ofF-diagonal long range order^^ to a mixture of two distinct types 
of bosonic atoms. Let 

.Na Nb 

7^a(x-x') = y / ]^dxi ]^dX/ *G(x,X2,...,XAf^;Xi,...,XAr3)*G(x',X2,...,XAf^;Xi,.. . ,Xnb) (A8a) 

•^ i=2 1=1 

and 

. Na Nb 

7^B(X - X') = y / n '^^^ n "^^i *g(xi, . . . , x^^; X, X2, . . . , Xjv^) *g(xi, . . . ,xw^; X', X2, . . . , Xjv^) . (A8b) 
•' i=i 1=2 

We can write these expressions in terms of the classical liquid correlations: 

Na Nb 

g-(7(x2,...,xiv^;Xi,...,Xiv^) 



7^A(x-x')= V—^ [f\dx^Y\dXi 

Zna.Nb J t^ fj-^ 



Xg-^[2:,'!4^AA(|x-X,|)+2:;;f,y^B(|x-X,,|)] ^-i[T.'^3,VAA{\^'-^,\) + T.J^^VAB{\^'-^j\)] .^g^ 

.Na Nb 



^ W^-UV, /TT^^^ TTdX,PG(x2,...,XAr,;Xi,...,X, 

Zna,Nb J V; 7_i 



^g-^[2:f=l!^AA(|x-X,|)+2:"fiVAB(|x-X,,|)] g-|[Ef=lVAA(|x'-X,|) + i:"fiVAB(|x'-X,,|)l^^QS 
= V ^^A-l,iVB ^g-l$A(x) g-^$A(x')^ (^^^-) 

where we defined the potential caused by the Na — 1 type A particles and the Nb type B ones as 

$a(x) - fdr [Vaa{\^ - r\)pA{r) + Vab{\^ - r|) psir)] , (A12) 

with the densities given by pA(r) = J2j=2 ^(''^ ^ ^j) ^^^ Pb(j) = J2j=i '^('^ ^ ^j)- 
Similarly, one obtains 

7eB(x-xO = y ^;^"-^""^ ^g-i*«(x) g-i$«(x')^ (A13) 

^AfA,JVB 

where the potential caused by the Na type ^ particles and the Nb — 1 type B is given by 

$b(X)= /"dr [y4B(|X-r|)pA(r) + VBB(|X-r|)pB(r)] . (A14) 

One can define the condensate fraction of type A bosons as ua = TIa{oo) = lim|x-x'|^oo T^a{^ — x'), and so 

„^ = y^^d^M^ (e-i*-(o))2, (A15) 

Zna^Nb 

where we assumed that the two-point correlation function factorizes, 

(e-^*^W e-^*^^''')) ^ (e-5*A(x)^ (e-^*^^''')) ^ (£-5*^(0)^2 ^ ^^^g) 

and also that the one-point function is translational invariant. Notice that the factorization assumption should be fine 
for short-ranged potentials or potentials that decay sufficiently fast, but there may be certain (possibly pathological) 
potentials for which it may fail. 

One can check that this definition of the condensate fraction of A bosons is such that ua < 1 as follows. First we 
use Holder's inequality to obtain 

(e-^*^ W e-^*-^(''')) < Y^(e-*A(x))(e-*A(x')) = (e-*^(o)) , (A17) 

from which we conclude that TZa{^ — x') < TZa{0) = 1, the last equality following trivially from the normalized 
wavefunctions (plus translational invariance). Thus, ha < 1 follows from the definitions above. Similar results apply 
to Ub- 

Notice that TZa{0) = 1 is a simple way to obtain (e^*-*^"') = y ^ "a.^b . .^^g .^^jjj ^gg ^j^jg equality below. 
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a. Fixing the chemical potential 

If we work at fixed chemical potential, we can relate the condensate fraction of bosons A and B to the fugacities 
zab a-iid the densities pa.b a-s follows. The grand canonical partition function for the binary mixture is given by 

Na-Nb 

At the saddle point, with particle numbers iV^ and Ng, one has the relations 

Zm* _i jv* = Zm' m' = Zm' m* _i . ( A19) 

"a '-'" B "A'"b „^ " A'" B ^ ^ ' 

ZA Zq 

Substitution of these relations in the equations for ua^b allows us to write 

UA^— (e-^*-^(°))2 (A20a) 

PA 



and 



nB = — (e-^*-(°V- (A20b) 

PB 

Notice that Eq. (A19) also allows us to write 

l^e-^Aio)^ ^ ]_ ^Na,n^ ^ PA ^^^ (e-^«(°)) = l ^^^'^" ^^ (A21) 

V Zna-i,Nb ^a V Znj^^n^^i zb 

at equilibrium. 

These relations allows us further to put lower bounds on the condensate fraction. To do so, consider without loss 
of generality potentials such that ^a,b{^) > 0; one can always do so by shifting the energies by a constant value such 
that the potentials Vaa, Vbb and Vab are non-negative (and using that the densities pa.b ^0). In this case, we can 
write 

(e-5*^.i5(o)) > ^g-*^,B(o)^ ^ PA^ (A22) 

ZA,B 

and thus 

nA.B > ^^^ . (A23) 

ZA,B 

For a hard sphere potentials, the equality (e~2*A,B(o)^ _ ^g-*A,B(o)^ jg satisfied. Thus, hard spheres have the 
lowest possible condensate fractions respecting the bound Eq. (A23). 
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crucial difference is that in Eq. (27) the Pa are normalized to f and the weights Wa add up to one, while this is not the case 
in Eq. (30). 

There is only one step that is not detailed in Refs. 33,34 so we sketch it here for completeness, and it is the proof that 
p(0) = - 4-|3/2 - In the replica formalism g?(x,x') — p~^'}2ii (^(^ ~ Xi')5(x' — Xj)), where a ^ h are two different 
replicas. In the language of Refs. 33,34, replicas are arranged in molecules labeled by the index i = 1, • • ■ , A''. Then \i i ^ j 
particles in different replicas are constrained to be at distance [x — x'| > a and do not contribute to p(x, x). The only term 
that contributes is the one for i = j, and from Eq. (4) in Ref. 34 we have p(x, x') = p~ X]i=i(^(^ ~ Xi)(5(x' — x^)) — 
p^^ J dx3 ■ ■ ■ dxm p(xi, ■ ■ ■ , Xm). Theu the result follows from Eqs.(ll) and (12) of Ref. 34. 



